arXiv:l 502.06325v 1 [nucl-th] 23 Feb 2015 


Cluster states and isoscalar monopole transitions of 24 Mg 

Y. Chiba and M. Kimura 

Department of Physics, Hokkaido University, Sapporo 060-0810, Japan 
(Dated: February 24, 2015) 

Background: Isoscalar monopole transition has been suggested as a key observable to search for exotic cluster 
states. Recently, the excited 0 + states with strong isoscalar monopole transition strengths are experimentally 
reported in 24 Mg, but their structures are unrevealed because of the lack of theoretical analysis. 

Purpose: Study structure of the excited 0 + states of 24 Mg populated by isoscalar monopole transition from the 
ground state and identify their cluster configurations. 

Method: The 0 + states of 24 Mg and their isoscalar monopole transition strengths from the ground state are 
calculated with antisymmetrized molecular dynamics combined with generator coordinate method using Gogny 
D1S interaction. 

Results: The calculated isoscalar monopole strength function shows reasonable agreement with experiment and 
is consistent with other theoretical calculation. The structure of the excited 0 + states with pronounced isoscalar 
monopole transitions are analyzed. It is found that the Oj, and 0^ states have mixed nature of mean-field 
and cluster, and that the 0^“ state is dominated by 12 C+ 12 C cluster configuration. In addition, it is predicted 
that 5ce-pentagon+<a states appear around 23 MeV. 

Conclusions: The excited 0 + states which appear as the prominent peaks in the calculated strength function 
are associated with 20 Ne+o, 12 C+ 12 C and 5o-pentagon+o cluster states. 


Introduction. — Clustering is a fundamental degree-of- 
freedom of nuclear excitation. According to Ikeda thresh¬ 
old rule [U, the appearance of various cluster states is 
expected near the cluster decay thresholds. Clustering 
of p shell nuclei has long been studied and well estab¬ 
lished 01 including the dilute gas-like (^-cluster state 
of 12 C(02") 0-0. On the other hand, in the mid sd-shell 
nuclei, the existence of cluster states is not well estab¬ 
lished, although many interesting phenomena can be ex¬ 
pected. For example, in the case of 24 Mg, a variety of 
exotic cluster states is expected: In addition to the ordi¬ 
nary a cluster state ( 20 Ne-j-ad, 12 C+ 12 C molecular states 
of astrophysical interest [^-[l8|. 16 0+2<a clustering E§- 
[22] and 6a condensation [7, Hj. |24| are theoretically dis- 
cussed. However, their high excitation energies make it 
difficult to identify them experimentally. 

In this decade, it is found that the isoscalar (IS) 
monopole transition strengths between the ground and 
excited cluster states are considerably enhanced, and 
hence, it can be a good probe for highly excited clus¬ 
ter states. The discussion was initiated by T. Kawa- 
bata [25j and Y. Kanada-En’yo [26| on the enhanced IS 
monopole transition of 11 B between the shell model like 
ground state and the 3/2^ state with pronounced 2 a + t 
cluster structure. T. Yamada et al [23, |28[ proved the 
mechanism of the enhanced IS monopole transition us¬ 
ing cluster-model wave function. The ingredient of the 
enhancement is the fact the ground state has “ duality na¬ 
ture ” of the mean-field and clustering f29, 30]. The dual¬ 
ity nature implies that the degrees-of-freedom of cluster 
excitation is embedded in the ground state even if it has 
a pure shell-model structure. It was shown that the IS 
monopole transition operator can activate this degrees- 
of-freedom of clustering. As a result, the excited clus¬ 
ter states can be strongly populated by the IS monopole 
transitions. In fact, the enhancements of IS monopole 


transition strengths are observed in p-shell nuclei such as 
n B [25[|, 12 C [3l[ and 16 0 [32[, and they nicely coincide 
with the cluster states predicted by theoretical calcula¬ 
tions. Thus, IS monopole transition is a promising probe 
for highly excited cluster states. 

Recently, the excited 0 + states with strong isoscalar 
monopole transition strengths are experimentally re¬ 
ported in 24 Mg [ 33 ], but their structures are ambigu¬ 
ous. Therefore, in this study, we aim to clarify the rela¬ 
tionship between those excited 0 + states and clustering. 
For this purpose, we calculate the excited 0 + states of 
24 Mg and investigate their clustering and IS monopole 
transition strengths from the ground state by the anti¬ 
symmetrized molecular dynamics (AMD), which has suc¬ 
cessfully described a variety of structure of p-sd-pf -shell 
nuclei [35l 137] including of the low-lying states of 24 Mg 
Q • To describe the cluster states and single-particle 
states including giant monopole resonance (GMR) simul¬ 
taneously, we introduce the constraint on the harmonic 
oscillator quanta and perform the generator coordinate 
method (GCM) with a large number of basis wave func¬ 
tions. 

Formalism. —We employ the microscopic Hamiltonian, 


# = Y ^ - tc.m. + F.VNNjij) + YFcfat (!) 

i i<j i<j 

where vn and vc stand for the center-of-mass ki¬ 

netic energy, Gogny DIS effective NN interaction [39j and 
the Coulomb interaction approximated by a sum of seven 
Gaussians, respectively. The AMD variational wave func¬ 
tion used in this study is an antisymmetrized product of 
the single particle wave packets projected to the positive- 
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parity state, 

* + ='^P L Mvx,V2,---'Pa}, (2) 



0 (ciiX't + hxi) ® (neutron or proton), (4) 

where the single-particle wave packet c Pi is represented 
by a deformed Gaussian wave packet J40|, and the vari¬ 
ational parameters z/ CT , ai and bi are determined by 
the energy variation. 

To deal with the low-lying quadrupole collective states 
and highly excited cluster states simultaneously, we in¬ 
troduce two different constraints in the energy variation. 
The first is imposed on the nuclear quadrupole deforma¬ 
tion parameters /3 and 7 to describe the low-lying collec¬ 
tive states, and we denote the set of the wave functions 
obtained with this constraint as . As the second con¬ 
straint, we extend the method used in Ref. j4l| and 
impose the constraint on the expectation values of the 
harmonic oscillator quanta N x , N y and N z , which are 
defined as the eigenvalues of the 3 by 3 matrix, 

A 

N ar = ($ + | X]4(*)«t(*)I $+ )) <T,T = X,y,Z. (5) 
i= 1 


various 0 + states from the low-lying to the highly excited 
states (GCM), 

*n + = E 9inP J= °$t, (8) 

where P J=0 is the projector to J = 0 state. We 
superposed 524 basis wave function 4>+ in total, and 
solved the Hill-Wheeler equation to obtain the eigenen- 
ergies E n and wave functions of the ground and 
excited 0+ states. To discuss the intrinsic structure 
of the 0 + states, we also calculate the overlaps be¬ 
tween each 0 + state and the basis wave functions, 

i« + \p j =°<f>+)\ 2 /(p j = o $i\p j = o $t). 

Using the wave functions of the ground and excited 
0 + states directly, we derived the IS monopole matrix el¬ 
ements M n {IS 0 ), reduced transition strengths B(ISO), 
strength function S(E X ) and the energy non-weighted 
and weighted sums rrik with k = 0,1, 3, 

M n (IS0) = {*° n + \d IS o\* 0 g + s.), (9) 

B(IS0-,g.s. -> 0+) = \M n (IS0)\ 2 , ( 10 ) 

S(E X ) = E \M n (ISQ)\ 2 E' n 6(E’ n - E x ), (11) 

n 

n OO 

m k = / dE x P\M n (IS0)\ 2 E^6(E' n -E x ), (12) 

^ 0 n 


Here a T (i) is an ordinary annihilation operator of the 
harmonic oscillator acting on the ith nucleon, and the 
oscillator parameter Tiuo is estimated from the ground 
state radius and set to 12.6 MeV. As a measure of 
the particle-hole excitation, we introduce the quantity 
AN = N x + N y + N z — TVo where Ab is the lowest 
Pauli-allowed value equal to 28. Under the condition 
of the AN = 0 , 2 ,4, 6 or 8 , we put the constraints for 
all possible even integer values of N x ,N y and N z . In 
other words, roughly speaking, we searched for the vari¬ 
ous many-particle-hole configurations within 8Tuj excita¬ 
tion. We denote thus-obtained set of the wave functions 
as $£ n . 

We further introduce an additional set of the basis 
wave functions defined as, 

*iso = (l - e-» 6lS0 ) $7 * fiOiso^, (6) 

A 

OlSO = E ( r * “ r c.m.) 2 , (7) 

i=1 

where p is arbitrary small real number, Oiso is the IS 
monopole operator and r c . m . is the center-of-mass coor¬ 
dinate. By definition, the set of the wave functions 
describes lplh (2fhj) excited states built on by the 
IS monopole operator. The similar method was also used 
in Ref. [42l. 

Those three sets of wave functions 4>^ , an d ®tso 
are projected to the J n = 0 + and superposed to describe 


where E' n stands for the excitation energy of the nth 0 + 
state, i.e. E' n = E n — E g . s .. 

Results of the energy variation .— Figure [T| (a) and 
(b) show the typical configurations obtained by the con¬ 
straint on the quadrupole deformation. After the GCM 
calculation, they become the dominant component of the 
ground and Oj states, respectively. The centroids of the 
Gaussian wave packets are gathered around the center- 
of-mass, describing triaxially deformed mean-field config¬ 
uration. As already discussed in our previous work |38j |. 
the constraint on the quadrupole deformation generates 



FIG. 1. Intrinsic density distributions at the z = 0 plane ob¬ 
tained by constraint on the matter quadrupole deformation 
parameters ((a) and (b)) and the expectation values of the 
harmonic oscillator quanta ((c)-(h)). The crosses in each fig¬ 
ure show the centroids of Gaussians describing nucleons. The 
contour lines are plotted in the interval of 0.02 fm -3 















3 


TABLE I. Calculated energy weighted sums m\ and ml in 
fraction of the EWSR and the centroid energies of GMR 
(ml/rriQ and yjm^/ml) in MeV, where rag, m* and are 
the sums between E x = 9 and 40 MeV excluding the 0^ state. 


basis set 

mi 

ml 

ml /mo 

\Jmltml 

(a) <L/3 7 

35 

26 

20.3 

24.2 

(b) $/3 7 +4 >/,so 

116 

101 

25.6 

29.3 

(c) 4> / 3 7 +4>/ ( 5 o+Tav 

103 

90 

22.2 

25.2 

exp. [45-47] 


82 ± 9 

21-9±g;l 

24.7 ±g;l 

QRPA [48] 


94 

20.57 



deformed mean-field configurations @ 0 ] , but no clus¬ 
ter configuration. 

The use of the constraint on the harmonic oscillator 
quanta generates various kind of cluster configurations as 
well as single-particle excited configurations with approx¬ 
imate ANTiuj excitations which lie energetically above 
the energy surface of and are not accessible by the 
constraint on the quadrupole deformation. The panels 
(c)-(h) show the examples of thus-obtained cluster wave 
functions, and they are the dominant component of the 
excited 0 + states corresponding to the prominent peaks 
in the IS monopole strength function S(E X ). By the con¬ 
straint of AN = 2 , 20 Ne + <a and 12 C + 12 C cluster states 
start to appear. As AN increases, the inter-cluster dis¬ 
tance grows and the orientation of cluster changes de¬ 
pending on the combination of N x , N y and N z . For exam¬ 
ple, the panels (c) and (d) show the 20 Ne + <a cluster con¬ 
figuration with AN = 4 and 8, which mainly contribute 
to the 0^ and O 5 " states, respectively. They have differ¬ 
ent orientation of 20 Ne cluster and inter-cluster distances 
(distance between the centroids of Gaussians describing 
clusters) are 3.0 and 3.3 fm, respectively. The panels (e), 
(f) and (g) show 12 C+ 12 C cluster states with AN = 6, 
6 and 8, respectively. They have different orientations 
of the oblately deformed 12 C clusters, and inter-cluster 
distances are 3.5, 3.5 and 4.0 fm. By further increase of 
A TV, very exotic cluster structure composed of 6 <a par¬ 
ticles appears. A typical example is shown in the panel 
(h) which was obtained by the constraint of AN = 8. In 
this configuration, centroids of Gaussians describing 5a 
clusters locate at the vertex of a pentagon with side of 1.5 
fm, and the last a cluster is 0.25 fm above it. After the 
GCM calculation, this kind of 5<a-pentagon + a configu¬ 
rations generate two 0 + states above 20 MeV. Thus, by 
increasing the number of particle-hole, 24 Mg is clustered 
as illustrated in Ikeda diagram [ 1 ]. 

IS monopole transition strengths. — The ground and 
excited 0 + states are calculated by the GCM with three 
different basis sets (a) <F/ 3 7 (b) <f>/ 3 7 + &iso ( c ) + 

&iso + n- The IS monopole transition strengths de¬ 
rived from these GCM wave functions are shown in Fig. 
[2j and their energy weighted sums and the centroid en¬ 
ergies of GMR are summarized in Tab. HI With only 
the basis set <1>£ 7 , the strength function (Fig. [2] (a)) 


fails to describe GMR, and the energy weighted sum mi 
amounts to only 35% of the sum rule (EWSR). Addi¬ 
tion of the basis set $>iso (Fig. [ 2 ] (b)) greatly improves 
mi value (116% of EWSR), but overestimates the ob¬ 
served GMR centroid energy j45l-l47t| because the GMR 
strength distributes widely in the region of E x > 30 MeV. 
The inclusion of the basis set &an yields the reasonable 
strength function as shown in Fig. [ 2 ] (c). Namely var¬ 
ious cluster and single-particle states with ANTiuo exci¬ 
tation described by &an lower the GMR position and 
enhance its strength. As a result, the strength function 
exhausts approximately 100 % of EWSR and plausibly 
agrees with the experimental energy weighted sum and 
the GMR centroid energy observed in the energy range of 
E x = 9 — 40 MeV. It is also noted that the quasi-particle 
random phase approximation (QRPA) with Gogny D1S 
interaction |39| also yielded similar values and qualita¬ 
tively agrees with our results and experiment with re¬ 
spect to the global structure of GMR. 

From the comparison between the strength functions 
shown in Fig. [2] (b) and (c), we also see that not only 
the GMR strength (E x > 18 MeV) but also the low- 
lying structure ( E x < 18 MeV) of the strength function 
is largely modified by the basis set &an- For example, 
note that the prominent peak at 15.3 MeV in Fig. [ 2 ] (c) is 
completely missing in Fig. [ 2 ] (b). Based on the analysis 
of the wave functions corresponding to those peaks, we 
conclude that several prominent peaks are attributed to 
the cluster configurations and suggest that the cluster 
states shown in Fig. |T] can be populated and observed 
by their enhanced IS monopole transition strengths. To 
see this point, we discuss the structure of the 0 + states 
relevant to the prominent peaks in S(E X ) in the following. 

Cluster states and their transition strengths .— The 
ground state is dominated by the mean-field structure 
and has the largest overlap (0.93) with the wave function 
shown in Fig. [1] (a). However, at the same time, it also 
has non-negligible overlaps with the cluster wave func¬ 
tions. It has 0.26 and 0.40 overlaps with 20 Ne+<a and 
12 C+ 12 C cluster states shown in Fig. [T| (c) and (e), re¬ 
spectively. This result means following two points. The 
first is that the cluster correlation exists even in the 
ground state. The binding energy of the ground state in¬ 
creases from 198.3 MeV to 199.2 MeV by including &an 
which indicates that the additional binding energy of 0.9 
MeV is brought about by the cluster correlation. Sec¬ 
ondly, it shows that the ground state has “ duality nature ” 
of the mean-field and clusters and that the degrees-of- 
freedom of cluster excitation are embedded in the ground 
state. This is an essential ingredient for the discussion 
of the IS monopole transition from the ground to the 
excited cluster states [27j. 

By including at, the low-lying mean-field states (the 
excited states having E x < 15 MeV in Fig. [ 2 ] (b)) are 
strongly mixed with 20 Ne + <a and 12 C+ 12 C cluster states 
and constitute the low-lying prominent peaks at 9.3, 11.7 
and 13.2 MeV in Fig. [ 2 ] (c), which correspond to the 
O 2 ", O 3 " and 65 “ states, respectively. In contrast to those 
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FIG. 2. The isoscalar monopole transition strength functions calculated with the basis sets of (a) <L/ 3 7 , (b) <f>/ 3 7 ± &iso and (c) 
4>^ 7 + &iso + 4>Aiv. The solid line in the right panel shows the strength function smeared by Lorentzian with 0.8 MeV width. 
The vertical dashed lines indicate cluster decay threshold energies which are located at the observed binding energies. 


mixed states, the 0g“ state at 15.3 MeV is dominated by 
the 12 C± 12 C cluster configurations. Furthermore, 5a- 
pentagon ± a cluster states configurations generates the 
0^5 and O 24 states at 21.8 and 24.9 MeV. 

The Oj state which appears as the lowest peak at 9.3 
MeV has the largest overlap (0.36) with the mean-field 
configuration of Fig. [1] (b) which has larger quadrupole 
deformation parameter /3 than the ground state. It can 
be regarded as the /3-band built on the ground band, 
and hence, has large IS monopole transition strength as 
listed in Tab. HD However, it also has 0.32 overlap with 
the 20 Ne ± a cluster configuration shown in Fig. [1] (c). 
Owing to this cluster correlation, it gains additional bind¬ 
ing energy of 1.8 MeV which reduces the excitation en¬ 
ergy from 10.2 MeV to 9.3 MeV. The 20 Ne±a cluster 
structure also constitute the Orj - state at 11.7 MeV by 
the mixing with mean-field structure. It has the largest 
overlap (0.30) with the configuration of the 20 Ne±a clus¬ 
ter configuration shown in Fig. [T] (d) and the comparable 
overlap (0.25) with the mean-filed wave function with 
(A 7 ) = (0.4,57°) in 

The 12 C± 12 C cluster configurations dominantly con¬ 
tribute to the Oj and 0^ states. The 0% state at 11.7 
MeV exhausts 3.7 % of EWSR and has the large over¬ 
laps with 12 C± 12 C, 20 Ne±a cluster and mean-field con¬ 
figurations. The overlaps are 0.21, 0.19 and 0.16 with 
12 C± 12 C and 20 Ne±a cluster and the mean-field con¬ 
figurations with (/?, 7 ) = (0.76,2.4°), respectively. In 
contrast to the above mentioned states, the 0g~ state 
at 15.3 MeV, which is close to the 12 C± 12 C cluster de¬ 
cay threshold energy, is governed by the 12 C± 12 C clus¬ 
ter configurations. It has 0.14 0.11 0.18 overlaps with 
the configurations of Fig. |T] (e), (f) and (g). The over¬ 
laps with other configurations are less than 0.09. The 
0 g" state has strong IS monopole transition strength and 
6.4 % of EWSR. It is noted that this state is completely 
missing in Fig. [2] (a), (b) and looks also missing in the 
QRPA calculation [48j , which is consistent with its strong 
12 C± 12 C cluster nature. Very interestingly, the 0 + state 


TABLE II. Properties of the 0 + states obtained by GCM cal¬ 
culation with 4>/3 7 , <f>/so and 4 >a n- E x , proton radius (r|), 
B(IS0) and E x B(IS0) are given in unit of MeV, fm, fm 4 and 
fraction of EWSR in percentag e, respe ctively. The values in 
bracket are the observed values [49l-[52| . The observed y/r% is 
deduced from the observed charge radius (5l[ and the proton 
charge radius [52]. 


State 

E x 

VW) 

B(IS0) 

E x B(IS0) 

Of 

0.0 

3.06 (2.93) 



Of 

9.3 (6.4) 

3.11 

122 (180 ± 20) 

6.1 (6.4 ±0.7) 

Of 

11.7 

3.08 

59.3 

3.7 

Of 

13.2 

3.06 

31.1 

2.2 

Of 

15.3 

3.11 

77.8 

6.4 

Ofs 

21.8 

3.14 

21.6 

2.5 

0f 4 

24.9 

3.28 

8.90 

1.2 


with enhanced IS monopole transition is observed at the 
12 C± 12 C cluster threshold energy in the 24 Mg(a,a') ex¬ 
periment [33j, which can be associated with the present 
0g“ state. We also note that the 0g“ state is the band- 
head of 12 C± 12 C cluster band which is a candidate of the 
observed 12 C± 12 C molecular resonances of astrophysical 
interest, for which we will discuss in detail in our forth¬ 
coming paper. 

Adding to those clusters, the 5a-pentagon±a clus¬ 
ter states appear as the 0^ 5 and 0 ^ states at 21.8 and 
24.9 MeV that exhaust 2.5 and 1.2 % of EWSR, re¬ 
spectively. The 0^5 state has the largest overlap which 
amounts to 0.43 with the 5a-pentagon±a cluster struc¬ 
ture shown in Fig. [T| (h). It also has non-negligible 
overlap ( 0 . 20 ) with the single-particle excited configu¬ 
ration described by &iso- The 0^ state has the largest 
overlap (0.26) with the similar configuration to that of 
Fig. [2(h) and has rather minor contributions from other 
configurations. One may attempt to associate these 5a- 
pentagon±a cluster states with dilute 6 a gas state anal- 
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ogous to the Hoyle state of 1 2 3 4 5 6 7 8 9 10 11 12 C(0^). Indeed, the recent 
Hatree-Fock-Bogoliubov calculation showed the possible 
existence of dilute not cluster states at very low density 
in N = Z nuclei 0 . However, we conclude that the 
calculated 0 f 5 and 0 j 4 states do not correspond to dilute 
a gas state suggested by the HFB calculation accord¬ 
ing to the following reasons. First, the radii of those 
states are not large enough to be a dilute a gas state 
and much smaller than those of Ref. [24j. Second, they 
are much more deeply bound compared to E x = 80 MeV 
reported in Ref. [24]]. Since AMD is free from the spu¬ 
rious center-of-mass kinetic energy and the parity and 
angular-momentum projections are correctly performed, 
the excitation energies of cluster states are greatly re¬ 
duced. In addition, the present and O 24 states are 
rather compact and have non-negligible interaction en¬ 
ergies between a clusters. Therefore, we conclude that 
those states have cluster structure but do not have di¬ 
lute gas nature. We conjecture that dilute 6a gas state 
will appear at higher excitation energy and to describe 
it, we will need to enlarge AN in order to include more 
spatially extended 6a configurations. Nevertheless, we 
emphasize that the exotic a cluster states are firstly ob¬ 
tained without a priori assumption on clustering in this 
study and it is shown that they are experimentally acces¬ 
sible via IS monopole transition from the ground state. 
Although they are embedded in the GMR energy region, 
we expect that they can be experimentally identified by 
their decay mode, because different from GMR, they will 
selectively decay through a particle emission. 

Conclusions. — In summary, we investigated the struc¬ 
ture of the excited 0 + states of 24 Mg and their IS 
monopole transition strengths based on AMD. The 
mean-field and cluster configurations of 24 Mg were ob¬ 
tained by the energy variation. In particular, by using 
the constraint on the harmonic oscillator quanta, the 


20 Ne+cq 12 C+ 12 C, and 5<a-pentagon+<a cluster config¬ 
urations were obtained without any a priori assumption 
on clustering. In addition, lplh (2fwj) excited configura¬ 
tions built by the IS monopole operator were also intro¬ 
duced as the basis wave functions of GCM. With these 
basis wave functions, the calculated 0 + states yielded rea¬ 
sonable IS monopole strength function. Namely, they 
exhausted almost 100 % of EWSR and reproduced the 
observed centroid energy of GMR. The result is also con¬ 
sistent with the QRPA calculation. 

We have shown that the several excited 0 + states with 
the enhanced IS monopole transitions are associated with 
20 Ne+cq 12 C+ 12 C and 5<a-pentagon+<a cluster configura¬ 
tions. The O 2 ", O 3 " and O 5 " states have the mixed nature of 
mean-field, 20 Ne+<a and 12 C+ 12 C cluster configurations, 
while the 0g“ state is governed by 12 C+ 12 C cluster config¬ 
uration. The 0g" state may be associated with the strong 
peak observed at the 12 C+ 12 C cluster threshold energy 
in the 24 Mg(cq a') experiment j33|. Furthermore, we pre¬ 
dicted that the 5<a-pentagon+<a cluster states appear in 
the GMR energy region. Even though they do not cor¬ 
respond to the dilute 6a gas state, it is emphasized that 
the exotic a cluster states were firstly obtained without 
any a priori assumption on clustering and were shown to 
to be experimentally accessible via IS monopole transi¬ 
tion from the ground state. We expect that the detailed 
comparison with the latest experimental data will reveal 
the exotic clustering of 24 Mg. 
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